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ABSTRACT 

At least two arguments suggest that the orbits of a large fraction of binary stars and extrasolar 
planets shrank by 1-2 orders of magnitude after formation: (i) the physical radius of a star shrinks 
by a large factor from birth to the main sequence, yet many main-sequence stars have companions 
orbiting only a few stellar radii away, and (ii) in current theories of planet formation, the region within 
~ 0.1 AU of a protostar is too hot and rarefied for a Jupiter-mass planet to form, yet many "hot 
Jupiters" are observed at such distances. We investigate orbital shrinkage by the combined effects 
of secular perturbations from a distant companion star (Kozai oscillations) and tidal friction. We 
integrate the relevant equations of motion to predict the distribution of orbital elements produced by 
this process. Binary stars with orbital periods of 0.1 to 10 days, with a median of ~ 2 d, are produced 
from binaries with much longer periods (10 d to ~ 10 5 d), consistent with observations indicating that 
most or all short-period binaries have distant companions (tertiaries) . We also make two new testable 
predictions: (1) For periods between 3 and 10 d, the distribution of the mutual inclination between 
the inner binary and the tertiary orbit should peak strongly near 40° and 140°. (2) Extrasolar planets 
whose host stars have a distant binary companion may also undergo this process, in which case the 
orbit of the resulting hot Jupiter will typically be misaligned with the equator of its host star. 

Subject headings: binaries: close — celestial mechanics — stars: planetary systems — methods: 
statistical 



1. INTRODUCTION 
1.1. Close binaries with tertiaries 

Close binary star systems (separation comparable to 
the stellar radii) are often accompanied by a third star. 
Such triple-star systems may be stable for long times if 
the system is hierarchical, that is, if the system consists 
of an "inner" binary (masses mi and 7712) in a nearly Ke- 
plerian orbit with semi- major axis »j„, and an "outer" 
binary in which 777,3 orbits the center of mass of the inner 
binary, with semi-major axis a out 3> »j„. If the forma- 
tion of the inner binary is independent of the formation 
of the outer binary, then the probability that a close bi- 
nary has a distant companion should be the same as the 
probability that an individual star has a companion in 
such an orbit. It is well known that the latter probabil- 
ity is substantial, about 2/3 for nearby systems whose 
brigh test component is a G-dwarf (|Duquennov fc Mayor! 
1991). Recent studies suggest that a similarly large 
fraction of contact binar i es have a third com ponent 
(jPribulla fc Rucinski|[200llD'Angelo et al.ll2006D . 

The data show, however, that characteristics of the 
inner and outer binary in hierarchical triples are not 
independent: the probability of having a third compo- 
nent turns out to be a func t ion of the period of the in- 
ner binary. iTokovinin et al.l (|2006f) found that 96% of a 
sample of spectroscopic binaries with periods less than 
3 d had a tertiary component, compared to only 34% 
of binaries with periods greater than 12 d. This sam- 
ple was not carefully selected to control biases, however 
such biases are unlikely to change the basic result, since 
the authors observe that all of the five binaries with 
period less than 10 d in the volume-limited sample of 
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iDucmennov fc Mayor! (| 1 9 9 If ) have at least one additional 
companion. 

A closely related observation is that the period dis- 
tribution of inner binaries in triple systems is quite 
differe nt from the period d is tribut ion of isolated bina- 
ries. ITokovinin fc Smekovl (|2002T) found a significant 
peak at about 3 d in the logarithmic period distribu- 
tion of inner binaries of triple systems, a peak that is 
not present in isolated binaries. Because of this fea- 
ture, the period distributions for binaries with and with- 
out a thi rd component are diffe rent with a significance 
of 0.999 (|Tokovinin et alj I2006T ). The period distribu- 
tion from time doma in surveys of eclipsing binaries also 
peaks at a few day s (jDevori 2005; P aczvhski et alJl2006t 
iDerekas et al.ll2007l ). but it remains unclear whether this 
may be attributed to a selection effect: binaries with 
large orbital periods have a lower probability of eclips- 
ing an d fewer eclipses per unit time which diminishes the 
signal (jGaudi et al.ll2005D . 

These observational results are surprising. In partic- 
ular, the median semi-major axis of the outer binary in 
systems with inner b inary period < 7 d (aj„ , < 0.07 AU) 
is a out ~ 70 AU ()Tokovinin et alj 12006 ) , similar to 
the m edian for all binaries in the iDuquennov fc Mavorl 
(|1991h sample. Why are the processes of star forma- 
tion correlated over a range of three orders of magnitude 
in scale? In this paper we explore the possibility that 
in some circumstances the distant companion enhances 
tidal interactions in the inner binary, causing its period 
to shrink to the currently observed value. 

1.2. Kozai cycles 

The study of the long-term behavior of three point 
masses interacting only through gravity has engaged 
physicists and mathematicians since the time of New- 
ton. In most cases, long-term stability requires that the 
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system is hierarchical (a ou t 3> a m ). An additional re- 
quirement for a stable hierarchical triple system is that 
the eccentricity e ou t of the outer binary cannot be too 
large, so that 7773 cannot make close approaches to mi 
or ni2. An equivalent statement is that the gravitational 
perturbations from 7713 on the inner binary must always 
be weak. However, even weak perturbations from the 
outer body can have important long-term effects on the 
inner binary. The simplest of these is precession of the or- 
bital plane, which occurs if the orbital planes of the inner 
and outer binary are not aligned. If the inner and outer 
binary orbits are circular, this precession is analogous to 
the precession of two rigid rings with the same mass, ra- 
dius, and angular momentum as the binary orbits: both 
the mutual inclination and the scalar angular momenta of 
the rings remain fixed, while the two angular-momentum 
vectors precess around the direction defined by the total 
angular-momentum vector of the triple system. 

An unexpecte d aspect of this behavior was discovered 
by iKozail (|1962f ). Suppose the inner binary orbit is ini- 
tially circular, with the initial mutual inclination between 
inner and outer binaries equal to imitiai- Kozai found 
that there is a critical angle i c such that if imitiai is be- 
tween i c and 180° — i c , then the orbit of the inner binary 
cannot remain circular as it precesses: both the eccen- 
tricity of the inner binary ei n and the mutual inclination i 
execute periodic oscillations known as Kozai cycles. The 
amplitude of the eccentricity and inclination oscillations 
is independent of the strength of the perturbation from 
the outer body, which depends on 7713 , a ou t , and e ou t , but 
the oscillation amplitude does depend on imitiai '■ for ini- 
tial circular orbits with imitiai = ic or 180° — i c , the max- 
imum eccentricity is 0, but if initial — 90° the maximum 
eccentricity is unity; i.e., the two inner bodies collide. 

Kozai cycles can be investigated analytically by av- 
eraging over the orbital phases of the inner and outer 
binaries ((Koiail H962T: IFord et al.ll2000h : this averaging, 
usually called the secular approximation, is justified be- 
cause the precession time is generally much longer than 
the orbital time of either binary. In the averaged prob- 
lem the semi-major axes of the inner and outer binary are 
both conserved. The analysis is simplest in the limiting 
case when a out 3> <Xi n (so that the perturbing potential of 
the outer body can be written in the quadrupole approx- 
imation) and the angular momentum of the outer binary 
is much greater than that of the inner binary (so that the 
orientation of the outer binary is a constant of the mo- 
tion). With these approximations, the following results 
hold, (i) The averaged quadrupole potential from the 
outer binary is axisymmetric relative to its orbital plane, 
(ii) The averaged problem can be described by a Hamil- 
tonian with one degree of freedom, (iii) The eccentricity 
of the outer binary is constant, (iv) The critical inclina- 
tion is i c — cos -1 \/3/5 — 39.2°. (v) If the inner orbit 
is initially circular the maximum eccentricity achieved 
in a Kozai cycle is e in , ma x = [1 - (5/3) cos 2 initial] 1 ^ 2 ■ 
(vi) Depending on the initial conditions, the argument of 
pericenter Ui n (the angle measured in the orbital plane 
between the pericenter of the inner binary and the or- 
bital plane of the outer binary) can either librate (os- 
cillate around 90° or 270°) or circulate. The system 
may remain at a fixed point with oJi n = 90° or 270° and 
e m = [l-(5/3)cos 2 7: / „] 1 / 2 if i c < i flx < 180° -i c . (vii) 



The only property of the Kozai oscillation that depends 
on the masses of the three bodies, their semi-major axes, 
or the eccentricity of the outer binary is t he period of the 
oscill ation, which is of order the timescale (Kiseleva et al. 
fl998h : 
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small-amplitude libration about the fixed point takes 
place with a period 
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(viii) Octupole and higher-order terms in the perturb- 
ing potential introduce a narrow chaotic zone around 
the separatrix between circulating and librating solu- 
tions as determined by the quadrupole approximation 
(jHolman et al.lll997f l 

Consider a sequence of triple systems in which the 
semi-major axis a ou t of the outer binary becomes larger 
and larger, while its mass, inclination, and eccentricity 
remain the same. The maximum eccentricity of the in- 
ner binary in the Kozai cycle will remain fixed, but the 
period of the Kozai cycle will grow as a^ ui .. This behav- 
ior will continue so long as the perturbation from the 
outer body is the dominant cause of apsidal precession 
in the inner binary orbit. Thus, weak perturbations from 
distant third bodies can induce large eccentricity and in- 
clination oscillations. However, small additional sources 
of apsidal precession in the inner binary — general relativ- 
ity, tides, the quadrupole moments of the two members of 
the inner binary, planetary companions, etc. — can com- 
pletely suppress Kozai oscillations caused by a distant 
third body if they dominate the apsidal precession. 

1.3. Tides 

The dissipative forces due to tides on the stars of the 
inner binary are only significant if the two stars are sepa- 
rated by less than a few stellar radii. Thus the tidal fric- 
tion on an inner binary with a semi-major axis of (say) 
0.25 AU, corresponding to an orbital period of 32 d if 
mi = T72 = l^f©, is normally negligible. However, if the 
amplitude of the eccentricity oscillation during a Kozai 
cycle is sufficiently large, the pericenter distance of the 
inner binary may become sufficiently small at some phase 
of the cycle that tidal friction drains energy from the or- 
bit, reducing the semi-major axis and thereby enhancing 
the friction, until the inner binary settles into a circular 
orbit wit h a semi-major axis of only a few stellar radii. 
Following lTokovinin et all (|2006l ). we refer to this process 
as Kozai cycles with tidal friction (KCTF). 

1.4. Previous work 

iHarringtonl (|1968f ) first suggested that KCTF might 
be an important evolut i onary mechanism for triple 
stars. iMazeh fe Shahaml (|1979f ) showed how the long- 
period perturbations of a third body could reduce 
the inner b inary's separa ti on on the tidal dissipation 
timescale. iKiseleva etall (|1998t) used KCTF to con- 
strain the strength of tidal dissipation in stars and fo- 
cused on the possibility that the inner binary of the Al- 
gol syst em might have shrunk significantly by this mech- 
anism. lEggleton fc Kiseleva -Eggleton (2001) presented 
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the differential equations that we use here to model 
KCTF, which parametrize the extra forces (tidal fric- 
tion, quadrupole from a distant companion, etc.) acting 
on a binary orbit. 

Kozai cycles have a l so be e n studied in the plan etary 
context. iMazeh et all l| 1997ft . Ilnnanen et~aT1 ((19971 ). and 
iHolman et al.l (|1997f ) all suggested that the large ec- 
centricities observed in many extrasolar planet orbits 
could be explained by Kozai cycles if the host star 
were a member of a binary system; however, this hy- 
pothesis leads to two predic tions that are not verified 
(jTremaine fc Za kamska 2004) : (i) high-eccentricity plan- 
ets should mostly be found in binary systems; (ii) multi- 
planet systems should have low eccentricities, since their 
mutu al apsidal preces si on su ppresses the Kozai cycle. 
Also, Ta keda fc Rasiol (|2005ft argued that even if ev- 
ery extrasolar planet host has an undetected compan- 
ion, Kozai cycles alone cannot explain the distribution 
of observed eccentricities. (They found that too many 
planets remained on nearly circular orbits, but that re- 
sult is sensitive to the initial eccentricity assumed since 
the time spent near the unstable fixed points at = 0, 
±45°, ±135° dominates the period of the Kozai cy- 



cle, and thi s time depends stron gl y on the initial eccen- 
tricity; see Ilnnanen et~all 119971 ) iKiseleva et alj ((19981 ) 
speculated that KCTF might explain the presence of 
massive planets on orbits close to their parent stars 
( "hot Ju piters"), a poss i bility that we re-examine in this 
paper. [Wu & Murray (2003]) elaborated the specula- 
tion of IKiseleva et al.l ((1998ft for the orbit of the planet 



HD 80606b, which is unusual because of its very large 
eccentricity (e = 0.93) and small pericenter distance 
( q(l - e) = 33 AU ) . 

iBlaes et al.l ((20021 ) have examined a process similar 
to KCTF for triple black-hole systems that might be 
found in the centers of galaxies; here the dissipative 
force is gravitational radiation rather than tidal friction 
but much of the formalism is the same. A major dif- 
ference is that tidal friction vanishes in a binary with 
a circular orbit and synchronously rotating stars, while 
gravitational radiation does not. Therefore the end-state 
of a black-hole triple subject to the analog of KCTF is 
a merger, leaving a binary black hole (which may be- 
come unbound by the gravitational radiation recoil of 
the merger). In the present application, however, gravi- 
tational radiation is negligible. Therefore the stellar and 
planetary cases offer opportunities to compare the orbits 
of the observed systems to the distributions predicted by 
the theory, which this paper quantifies. 

1.5. The plan of this paper 

In this paper we shall model KCTF using (i) the secular 
approximation for orbital evolution; (ii) the quadrupole 
approximation for the tidal field from the outer body; 
(iii) a simple model for apsidal precession that includes 
the dominant general relativistic precession term and the 
quadrupole distortion of the stars of the inner binary due 
to tides and rotation; (iv) a simple model for orbital de- 
cay due to tidal friction; (v) the assumption that the 
outer binary contains most of the angular momentum in 
the system. Our main goal is to characterize the statis- 
tical properties of the binary systems that result from 
KCTF. In <2] we describe the equations of motion we 
use to evolve hierarchical systems under KCTF. In Sj3]we 



derive a Hamiltonian formalism and the conserved quan- 
tities of the system and thereby offer analytical insight 
into the behavior of Kozai cycles. We evaluate the ef- 
fect of tidal friction on a population of isolated binary 
stars in Sj4] We numerically integrate a large number of 
triple systems in fj5j compare the final conditions to an 
empirical period distribution, and give a prediction for 
the mutual inclination distribution. In Sj6]we verify that 
the mutual inclination distribution, once established by 
KCTF, will persist for the lifetime of the system, at least 
for isolated triple systems. We apply our results to hot 
Jupiters orbiting one component of wide binary systems 
in SjT] We discuss applications and possible extensions of 
our analysis in £|S]and restate our main conclusions in 

2. EQUATIONS OF MOTION 

The differential equations that govern the inner bi- 
nary's orbital parameters and the two stars' spin param- 
eters were presented by lEggleton fe Kiseleva-Eggletonl 
((2001ft and are recalled here. The equations take into 
account stellar distortion due to t ides and rotation, tida l 
dissipation based on the theory of Egglet on et al.l |l998), 
relativistic precession, and the secular perturbations of 
a third body (averaged over the inner and outer Ke- 
plerian orbits). The orientation of the inner orbit is 
specified by its Laplace-Runge-Lenz vector, ei„, whose 
magnitude is the inner eccentricity ei n and whose di- 
rection is towards the pericenter of the orbit of m 2 
about mi. Also needed is the specific angular momen- 
tum vector hi n of the inner orbit, whose magnitude is 
[G(mi + TO2)a.;„(l — eL)] 1 ^ 2 , where G is the gravita- 
tional constant. The vector q„ = h,;„ x e,- n completes 
the right-hand triad of unit vectors (e^, qi ra , hj n ). Each 
of the stars of the inner binary also has a spin vector ttj . 
J = 1,2. 

The evolution equations are: 



1 de m 

3 in, dt 



= {Z 1 +Z 2 +Z GR )q in -(Y 1 +Y 2 )h in -(V 1 +V 2 )&in 



(1 e iri) ['iSeq&in (45* ee Sqq)£[i n + Sgfrhin] , (3) 
1 dill. 



u = (Y 1 + Y 2 )e m -(X 1 +X 2 )q m -(W 1 + W 2 )h l . 
kin dt 

+ (1 - e 2 m )S q he in - (4e?„ + l)S e hAin + 5e- n S e qh. in , 



dfii 

h— tt =^hin{-Y 1 e in + A"iq m + Wihf, 
at 



d(l 2 

'~dt 



">h in (-Y 2 e m + X 2 q m + W 2 hi, 



(4) 
(5) 
(G) 



Subscripts 1 and 2 refer to the masses mi and m 2 . The 
gravitational influence of the third body is described by 
the tensor S. We specify the orientation of the orbit of 
the third body using the triad (e out , q out , h out ) , defined 
by analogy with the inner binary's coordinate system, 
treating the inner binary as a point mass mi + m 2 at its 
center of mass. We have: 



C{8. 
C 



in a 3/i 

2tt 



3t 



(1 



out ,m h out ,n) 
2 s-1/2 



(7) 
(8) 
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Here the directions m,n£ (ei n ,qi n ,hi n ) are along the 
basis vectors of the inner orbit, and h out . m and h out . n are 
direction cosines of the outer orbit's angular momentum 
along the inner orbit's coordinate directions. Pj„ and 
Pout are the periods of the inner and outer orbits, respec- 
tively. V and W are dissipation rates for e and h, respec- 
tively. The vector (X, Y, Z) in the (e in-, Qm : hm) frame is 
the angular velocity of that frame relative to the inertial 
frame. Zqb, is the first post-Newtonian effect of relativ- 
ity which causes the pericenter to precess. The stars of 
the inner binary have moment of inertia I\ and I2, and 
fi = ui\mil {rn\ + m-i) is the reduced mass of the inner 
binary. The functional f orms of the V, W, X, Y, and 
Z term s were modeled in lEggleton fc Kiseleva-Eggletonl 
(|2001h and for convenience are repeated in the appendix. 

To check the code we use to int egrate these equations , 
we have replicated the results of IWu fc Murravl (2003), 
who studied KCTF for the planet HD80606b, whose host 
star has a binary companion with a ou t ~ 10 3 AU. We 
recalculate their Figure 1, using their specified param- 
eters and initial conditions (Fig. [1]). We started both 
the planet's and the star's obliquity (the angle between 
hin and fii or O2) at 0°. Prominent eccentricity oscil- 
lations are seen in panel (a). The energy in the planet's 
spin was transferred to the orbit, increasing the semi- 
major axis for the first 0.1 Gyr (panel b). A rotational 
equilibrium was reached which matched the spin angular 
velocity to that of the instantaneous orbit angular ve- 
locity when tides were the strongest, i.e., for pericenter 
passage during the high-eccentricity part of Kozai cy- 
cles. As dissipation shrunk the semi-major axis, pericen- 
ter precession became gradually dominated by relativity 
rather than by the third body. Thus, although the orbit 
initially had a librating pericenter, it started circulating 
as the eccentricity passed close to zero at 0.7 Gyr (panel 
a). During the Kozai cycles, the orbital angular momen- 
tum of the planet is generally quite misaligned from the 
spin angular momentum of the star (panel e); i.e., the 
stellar obliquity tf> is substantial. After ai n equilibrates, 
a large amount of misalignment will generally remain for 
the main-sequence lifetime of the star, although slight 
movement towards alignment occurs because of dissipa- 
tive tides being raised on the star by the planet. The 
tidal field of the planet applies a torque to the rotational 
bulge of its host star, which causes the spin axis to pre- 
cess; compared to the fixed reference frame of the third 
body's orbit, the host star's spin axis oscillates with large 
amplitude (panel f). To conserve angular momentum, 
the planet's orbit likewise precesses, but on a consid- 
erably smaller scale — it is the few degree wiggle in the 
inclination after migration (after 2.8 Gyr in panel c). We 
will consider these alignment issues in <JB] and i)7.4l 



3. ANALYTIC UNDERSTANDING OF KCTF 

In this section we seek an approximate analytical the- 
ory of KCTF. To lowest order in ai n /a out , the ter- 
tiary component presents a quadrupole tide to the in- 
ner binary. As a,i n shrinks by KCTF, this approxima- 
tion will become better and better. The instantaneous 



quadrupole-order Hamiltonian is (|Ford et alJl2000j) : 

Gra\m,2 G{m\ + 1712)1713 



2&i n 

Gmim 2 TO3 rf n 
mi + m 2 rl m 



2a 



P 2 (cos$) 



■Fa 



(9) 
(10) 



where r^ n is the vector from mi to m%, r ou t is the vector 
from the inner binary center of mass to 7713, 4> is the angle 
between rj„ and r out , and = \x 2 — |. 

This quadrupole-order potential may be integrated 
over the unperturbed motion in both orbits to remove 
short period terms, which depe nd on the two orbits ' 
mean anomalies. However, as in llnnanen et al.l (|1997h . 
we shall work in the approximation that the outer or- 
bit contains essentially all the angular momentum and 
hence is fixed, therefore the mutual inclination i is the 
same as the inclination relative to the total angular mo 
me ntum (ij n = i. i nui ■ = 0). The averag ed Hamiltonian 
is (|innanen et al.lll997t iFord et aLlfeOOOh : 



(7) 



Gm\m2 G(mi + 1122)1713 



2a in 



Gmim27ri3 



2a a 



+ (F g ) (11) 



m!+m 2 8a 3 out (l - e 2 out ) 3 / 2 



x (2 + 3e 2 „ - (3 + 12e?„ - 15e?„ cos 2 u in ) sin 2 tj , (12) 

where the orbital elements are referenced to the plane of 
the outer binary. The canonical Delaunay variables for 
the inner binary are the mean anomaly U n , argument of 
pericenter Wi n , and longitude of the ascending node flm , 
along with their respective canonical momenta Li n — 
m 1 m2^Ga in /(iTii + m 2 ), G in = L ln ^J\~ e\ n = (j,h in , 
and Hi n = Gi n cosi. The canonical equations of motion 
for the inner orbit are: 



a, 



d{F) 

dLi n 


Lin — 


d{T) 

dlin 


d(T) 

dGi n 


Gin = 


du> in 


dH in 


Hin — 


d(T) 
BO, 



(13) 
(14) 
(15) 



with corresponding equations for the evolution of the co- 
ordinates and momenta of the outer variables. 

The averaging procedure removes the Hamiltonian's 
dependence on the mean anomalies /,•„ and l ut, so Ly n 
and Lout and thus the semi-major axes of the inner and 
outer orbits are conserved. The average Hamiltonian (T) 
is conserved, because the Hamiltonian is independent of 
time; moreover (J- q ) depends only on and the con- 
served semi- major axes (eq. [H]), so {T q ) is also con- 
served. Another conserved quantity is Hi n , because the 
Hamiltonian is independent of its canonical conjugate, 
Qi n . A final conserved quantity is the eccentricity of the 
outer binary, e out , because the Hamiltonian is indepen- 
dent of u> ou t ■ 

We may write dimensionless versions of the conserved 
quantities as: 
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Fig. 1. — The evolution of a planet initially in an orbit with = 5 AU, i = 85.6°, ei n = 0.1, uii n = 45°, as a hypothetical progenitor 
to HD 80606b. The spins of both the planet and its host star were initialized with zero obliquity. The stellar companion was assumed to 
have 7 713 = 1.1m©, a ut = 1000 AU, and e out = 0.5. The values of structural parameters were the same as those used by Wu & Murray 
i2003h . and the viscous times were tv iS tar = 50 yr and ty p [ anBt = 0.001 yr (see appendix). Energy dissipation is dominated by the planet. 
The diamonds mark the current position of HD 80606 along this possible evolution. 



+ (3 + 12e 2 „ — 15e 2 „ cos 2 u m ) sin 2 i (16) 
ff'^l-e^oosi. (17) 

iKiseleva et all (|1998l) have given a different conserved 
quantity which is simply a combination of the above: 
{5-3H' 2 +F')/3. 

These constants can be evaluated with the initial set 
of elements {ei n ,i n iUal, Uin,*nitial, iiniUal)- It is then pos- 
sible to compute analytically the values of e; n and i at 
any value of u>i n accessible to the system. Let us restrict 
our discussion to iinitiai > 90° (for retrograde systems, 
identical behavior of the inclination, mirrored across 90° , 
results). The system attains a maximum ej„ and mini- 
mum i when u>i n — 90° or 270°. Therefore, we have: 



e in>max = ({[(10 + 12H' 2 - F') 2 - 54CLff' 2 ] 1/2 

+8-12 J ff' 2 + F'}/18) 1/2 (18) 

W = co S - 1 [fl"(l-e 2 „ imax )- 1 / 2 ]. (19) 

If F' is initialized with = 0, the value of ei n _ max 
given by equation (|18p matches the value quoted in the 
introduction. 

Figure [5^ is a contour plot of the minimum mutual 
inclination attained as a function of initial eccentricity 
and argument of pericenter. For an initial eccentricity 
Gin,initiai < 0.2, the minimum inclination is within 2° 
of the critical inclination, i c = 39.2°. Integrating over 
a uniform initial distribution of arguments of pericen- 
ter and various initial eccentricity distributions gives the 
probability distributions of minimum inclination shown 
in Figure [2Jd. These results are only weakly dependent 
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on the initial inclination. 

We now argue that the minimum inclinations shown in 
Figure [5^ are very nearly equal to the final inclinations 
produced by KCTF. If we approximate the effect of tidal 
dissipation as acting only when the inclination is at its 
minimum (i.e., eccentricity maximum), then i m in is con- 
served between Kozai cycles (Fig. even though the 
constants F' and H' are not conserved in the presence 
of tidal friction. Eventually, all eccentricity is damped 
(after 4 Gyr in Fig. and the system finishes with an 
inclination nearly equal to its minimum inclination from 
the first cycle (Fig. [Tfc). The distributions of minimum 
inclination shown in Figure [2)3 will correspond to the fi- 
nal inclinations after KCTF if the above approximations 
are correct. Plotted for comparison is the inclination 
distribution of an isotropic distribution of triples, i.e., 
triples with inner and outer angular momentum direc- 
tions uncorrelated. The distributions clearly distinguish 
a population of triples whose inner binaries have shrunk 
by KCTF from a population with its primordial inclina- 
tions, if these are uncorrelated. 

This calculation will only be a faithful representation of 
the final inclination distribution if pericenter precession 
due to other causes is negligible compared to the pre- 
cession caused by the tertiary component. Under these 
conditions, the final distribution of mutual inclination 
shows strong peaks close to the critical angle for Kozai 
cycles: 39.2° or 140.8°, as seen in Figure [2b- For a more 
realistic calculation considering a population of triples 
and including the extra precession forces, see $5] and Fig- 
ure [71 

3.1. Kozai cycles in the presence of extra forces 

Because Kozai cycles are driven by the interplay be- 
tween the weak tidal forces from the outer binary and 
the shape of the orbit of the inner binary, they can 
be easily suppressed by other weak effects leading to 
pericenter precession in the inner binary. Mathemati- 
cally, the rate of change of eccentricity is proportional to 



(cleg) 



3Ch. 



out^e^out^q 



cx sin2wi„ (according to eq. [3] 
and eq. [7])- If is changing too fast, S eq will average 
to zero before ei n has time to grow. We consider four 
causes of extra pericenter precession — relativity, tides, 
stellar rotational distortion, and extra bodies. Here we 
characterize the first three effects with additional Hamil- 
tonian terms and calculate the associated apsidal motion 
u>i n . For extra Hamiltonian terms (Fextra), equation (I14| 
gives: 



^extra 



J 1/2 d(T extra ) 



dei 



(20) 



First, let us consider relativ ity. As in 

lEggleton fc Kiseleva-Eggletonl (|2001h . we only con- 
sider pericenter precession, which is the largest effect. 
Take an extra Hamiltonian term of 



(^Gfl}(etn) 



3G 2 ra\ra 2 (ra\ + 7712) 



1 



»2 \l/2 ' 
i) 



a m c2 (1 - ej 

(21) 

which depends on both L in and Gi n . The former de- 
pendence gives an additional contribution to the mean 
motion, which is incorrect but does not affect the secular 




Fig. 2. — (a) The minimum inclination during a Kozai cycle, 
starting with initial = 86° at various eccentricities e; n initial and 
arguments of pericenter ui in initial for the inner orbit (solid lines). 
These results depend very weakly on the initial inclination; for 
^initial = 70°, the contours for i min = 25°, 38°, 41°, and 65° (dot- 
ted lines) are only slightly different, (b) The inclination minimum 
distribution (probability density function), assuming uj in initial is 
uniformly distributed in angle. These distributions result from 
three different assumptions for the initial eccentricity distribution: 
circular orbits, uniform distribution in eccentricity e in initial? an d 
uniform distribution in e 2 . . , , (constant phase-space density on 

the energy surface). They were computed with initial = 86°, but 
all initial inclinations near 90° (i.e., all systems that evolve sub- 
stantially by KCTF) produce similar distributions. If KCTF seals 
in this inclination minimum (see text), the distribution of triples 
will show spikes at the critical angles i c ~ 40° and 180° — i c — 140° 
and very few systems near 90° . An isotropic distribution, with no 
correlation between the directions of inner and outer orbital an- 
gular momenta, is plotted for comparison. See also Figure [7p for 
the inclination distribution that results from integrating the full 
equations of motion for a population of triples. 



results. The latter gives 

3G 3 / 2 (m! 



^GR 



m 2 



)3/2 



4{ 2 c 2 (l - e?J 



(22) 



which is the standard expression for the rate of pericenter 
precession due to relativistic effects. 

Next let us consider how the non-dissipative tidal bulge 
contributes to the apsidal motion. At any moment, con- 
struct a spherical polar coordinate system centered on 
rn,2 with radius r 2 and a polar angle 9' measured from 
the vector fi n . The tidal potential of mi to lowest order 
in r 2 /r in is 



= -GmiP 2 (cos9')rli 



(23) 



The surface of m 2 is an equipotential. The corresponding 
distortion produces an external potential of 



52, ep 



= -2Gm Y P 2 (cos 9')r. 



-3-3 d5 



R 2 k 2 , 



(24) 



where Rj is the radius and fc , is the classical apsidal 
motion constant (|Russelllll928h of the jth body. Typical 
values of kj are 0.014 for stars and 0.25 for gas giant 
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planets. Back at mi (for 7-2 
force per unit mass is 

d(j>2, ep 



6' = 7i"), the extra 



dr 2 



-6Gmir iTI R\k 2 , 



(25) 



which can be integrated to find the potential associated 
with the tidal distortion of m 2 by mi'. 

0i,2 - -Gm ir ^Rlk 2 . (26) 
Notice this is a factor of 2 smaller than equation (I24| 
evaluated at the location of m\ because the derivative of 
02, ep is taken with respect to the spatial coordinate r 2 , 
but 01,2 is assembled by integrating: equation (|25|) with 
re spect to r in (see, e.g.. ISterndll939l and the Appendix 
of lEggleton fc Kiseleva-Eggletonll200l[) . By analogy with 
the theory of image charges in electrostatics, the effective 
potential is half the physical potential (P. Eggleton, pri- 
vate communication). From this potential we may form 
the instantaneous Hamiltonian Tride,2 = mi<Pi,2- After 
averaging over one orbit of the inner binary and account- 
ing for both stars, we have 

/T w v G 8 + 24eL + 3ei 

{■r Tide) {Bin ) 1 



Hn (1 

m\kiRl ■ 



- e 2 ) 9 / 2 
m\k 2 R\ 



(27) 



After converting to canonical variables, the equations of 
motion yield: 



^Tide : 



15(G(mi + m 2 )fl 2 8 + 12e 2 „ + e| 



13/2 



(1 



m 2 



mi 



— kiR\ + —k 2 R 



(28) 



Lmi m 2 
This expression is always positive, so tidal bulges always 
tend to promote pericenter precession and therefore sup- 
press Kozai oscillations. 

An analogous procedure gives the extra piece of the 
Hamiltonian arising from the rotational bulges of the 
stars of the inner binary. The instantaneous Hamilto- 
nian resulting from the quadrupolc field of the rotation- 
ally oblate 77i2 is: 



3~ Rotate, 1 ! — \k 2 



mi^&Rl 



P 2 (cos( 



(29) 



where 9 is the angle measured from the spin axis of 
m 2 . Averaging the result, accounting for both stars and 
putting Sli and il 2 in component form, 

77117712 



oi ai e }(ei n , Ui,Sl 2 ) 



6<4(1 



)3/2 



J=l,2 



(30) 



To use the canonical equations of motion, each f2j com- 
ponent must be converted to components in the inertial 
frame. After taking a derivative with respect to Gi n and 
converting back to the orbit frame, we have: 

(mi + m 2 ) 1 / 2 

^Rotate 



2GV2 a J/ 2 (l-4j; 



E 

3 = 1,2 



k R] 



{2Q 2 



n 2 ) 



2Qjh cot i(flj e sine 



ttjq COSU) in ) 



(31) 



Small additional bodies in the system, such as plan- 
ets, also affect the dynamics through their contribution 
to the apsidal precession rates of the stars in the in- 
ner binary. In binary stellar systems, planetary orbits 
may stay close to one of the stars (S-type orbits, for 
"satellite") or encompass the whole binary (P-type or- 
bits, for "planetary"). Stable S-ty pe or P-type orbits 
must obey certain stability crite r ia dHolman fc Wieg:ertj 
[19991: IMardling fc Aarsethl l200il: IMudrvk fc WuT l2006'l: 
crudely, these require that the semi-major axis a <ti ai n 
for stable S-type orbits, while a 3> a.i n for stable P-type 
orbits. In triple stellar systems, an analogous classifica- 
tion scheme can be worked out. Three different types of 
orbits may be stable: (1) S-type about any of the stars, 
(2) P-type with respect to the inner binary, but S-type 
with respect to the outer binary, and (3) P-type with 
respect to the outer binary. In the first case, the time- 
average of an S-type orbit in the equatorial plane of its 
host star will qualitatively act as an additional contri- 
bution to stellar oblatcness. Although large eccentricity 
oscillations of the inner binary would tend to destabilize 
S-type planets, the extra pericenter precession caused by 
the planet may suppress those oscillations: even a tiny 
planet may thereby be responsible for its own survival. 
In case (2) the outer binary can induce Kozai oscillations 
in the planetary orbit, while in case (3) the planet can 
induce Kozai oscillations in the outer binary. The con- 
tribution of such additional bodies to Kozai oscillations 
has been studied in the context of multi-planet S-type 
systems in binaries (llnnanen et ail 119971 : IWu fc Murray! 
120031 : iMalmberg et al.l [20061 ). but the huge parameter 
space of the general secular four-body problem has not 
yet been systematically explored. 

3.2. Modified eccentricity maximum 

When the additional forces described in §3.11 are 
present, we may compute e max as before, although the 
conserved Hamiltonian is now 

{Ftot) = {Fa) + {Fextra), (32) 

where {T extra ) = (Pgr) + (J^Tide) + (^Rotate)- If 
(^Rotate) contributes significantly, one must take into ac- 
count that each component of the spins (Flu, &2i) will 
be a function of time. 

Let us consider the case in which (F e xtra) is dominated 
by {Pgr) an d assume that the orbit is initially circular 
to evaluate (J-'tot)- Eccentricity maximum still occurs 
at bJi n = 90° or 270°, which simplifies the expression 
for (T q ) at eccentricity maximum. The conservation of 
{Ftot) and H' (eq. |17| ) gives an implicit equation for 



COS iinitial — 



2 

i7i, max ) 



r i/2 + ir 



tujgr\ 



(33) 



where the strength of relativity is parametrized by the 
product of the Kozai timescale, r (eq. [1]) and u>gr eval- 
uated at e.i n — (eq. [22]). We plot this function in 
Figure There is a maximum eccentricity that can be 
reached for a given amount of GR precession; beginning 
with mutual inclination of 90° and negligible eccentricity: 



= {l-[(* + §™ GJl | e . =0 ) 1/2 -l/2] 2 } 1 / 2 . (34) 
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Fig. 3. — The maximum eccentricity attained by systems with 
initially circular orbits of va rying initial inclinations, including rel- 
ativistic precession (eq. [33]). The curves are parameterized by the 
relative strength of relativistic precession to that of the tidal field 
of the third body: tu>gr\ c . _ , where r is defined by equation ifl) 

and u>gr (eq. |22| ) is evaluated at e; n = 0. As the timescale of rel- 
ativistic precession becomes as short as the timescale of precession 
induced by the companion, the eccentricity cycles are reduced in 
amplitude and the critical inclination — the largest iinitial f° r which 
the orbit remains circular — grows (eq. 1351 ). In contrast to the sit- 
uation without relativistic precession, a system in which the inner 
and outer binary orbits are initially perpendicular (iinitial = 90°) 
will not reach a radial orbit but attain a more moderate eccen- 
tricity (see eq. 1341 ). No Kozai oscillations occur for any initial 
inclination if tloqr\^ _ > 3. 



The critical inclination for eccentricity oscillations is 
also increased by relativity: 

cos 2 i c = 3/5 - (l/5)rw Gi? | e!ji=0 . (35) 

In summary, some hierarchical triples, even if they begin 
with perpendicular orbits, may avoid close encounters 
because relativistic precession suppresses K ozai cycles. 
A sim ilar analysis has been performed by iBlaes et al.l 
([20021 ). 

4. TIDAL SHRINKAGE IN ISOLATED BINARIES 

In this section, we select a sample of isolated binaries 
via Monte Carlo methods, then we follow how their or- 
bits evolve by tidal dissipation. This calculation provides 
a control sample for the numerically integrations of the 
following section, in which we follow the secular evolution 
of triple stars by KCTF. 

The initial orbital distr ibutions were taken f r om th e 
observed distributions of iDuquennov fe Mavorl (|1991f l. 
For simplicity, toi was set to M Q . Next, 7772 was 
chosen by selecting the mass ratio qi n — m 2 /77ii 
from a Gaussian distribution with mean 0.23 and stan- 
dard deviation 0.42, as found by IDuquennov fc Mavorl 
(|1991|) (negative results were resampled). The stars 
were given radii consistent with their masses for stars 
on the main s equence: R\ = and R2 — 

i? (m 2 /M ) - 8 (jKippenhahn fc Weigertj 11994ft . Their 
structural constants were set t o typ ical values given in 
lEggleton fc Kiseleva-Eggletonl (|2001h : ki = k 2 = 0.014 
(apsidal motion constants), tyi = tv2 = 5 yr (vis- 
cous timescales; see appendix), 1\ = 0.08Afi_R 2 and 
I 2 = 0.08M 2 i?2 (moments of inertia) . The initial spin pe- 
riods were each 10 d, with spin angular momenta aligned 
with the binary's orbital angular momentum. 



The orbital peri od was picked from the log- 
normal distribution of lDuquennov fc Mavorl (|1991l ) (with 
(logP[d]) = 4.8 and <7 logP [ d ] = 2.8). It is likely that a 
lower cutoff needs to be imposed because the large radii 
of pre-main-sequence stars likely preclude periods below 
a few days, and a cutoff of 0.2 d is seen empir ically for 
main-sequence binaries (jPaczyhski et al.l 12006) . We do 
not impose any cutoff at this stage, but we will explore 
such cutoffs when presenting the results. 

The distribution of initial eccentricities was a func- 
tion of period, in acco rd with observed binaries 
()Duquennov fc Mavorlll991h . For P < 1000 d, e was cho- 
sen from a Rayleigh distribution (dp oc e exp(— Ae 2 )de) 
with (e 2 ) 1 ' 2 = (2A)- 1 / 2 = 0.33; and for P > 1000 d, 
e was chosen from an Ambartsumian distribution (dp = 
2ede), corresponding to a uniform distribution on the 
energy surface in phase space. One might suppose that 
the observed eccentricity distribution is peaked towards 
small values for P < 1000 d solely because of tidal dissi- 
pation, but this cannot be the case. In particular, only 
one out of 12 isolated binaries wit h 11 d < P < 1000 d 
in the volume-limited sample of IDuquennov fc Mavorl 
(|1991f) has e > 0.5. However, tidal dissipation would 
only be able to affect the orbit of a main-sequence bi- 
nary with P 100 d if e > 0.8. Therefore, the small 
eccentricities of binaries with P < 1000 d cannot mostly 
result from the circularization of binaries with larger ec- 
centricity. 

A total of 7 x 10 4 such binaries were selected, and the 
evolution of P, e, fii and ^2 were followed using the 
equations of <J2] The integrations proceed very quickly, 
since the spins were assumed to be aligned with the 
orbit — only the magnitudes of the vectors in equations 
©-JS]) needed to be followed, not their orientation. The 
integrations were stopped when e < 10~ 3 or at 10 Gyr 
(roughly the main sequence lifetime of a l-M© star), 
whichever came first. 

In Figure [4] we plot the initial and final period distri- 
butions. The Gaussian tail at long periods is not plot- 
ted. Histograms corresponding to various initial period 
ranges are plotted in different shades. In particular, if 
the primordial period distribution were cut off below 6 d, 
both the initial and final histograms would consist only of 
the two darker shades. The unshaded region is probably 
not physical — main-sequence binaries with such periods 
would merge. 

A firm result is that tidal dissipation alone causes very 
little change to the period distribution of binaries. We 
turn now to KCTF. 

5. POPULATION OF TRIPLES 

In this section we calculate the semi-major axis and in- 
clination distributions that result from KCTF. We gen- 
erated systems by Monte Carlo methods, then integrated 
their equations of motion (Sj2J). 

5.1. Initial conditions 

The initial orbital distributions for both the inner and 
o uter binaries were taken from the observed distributions 
of lDuquennov fc Mavorl (|199lD . Both the orbital distri- 
butions of the binaries and the physical parameters of 
the stars of the inner binary are chosen as above (Q. 

For the tertiary, 7773 was determined by choosing q out = 
7713/(7711 + 7772) from a Gaussian distribution as for qi n . 
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P [d] 

Fig. 4. — The periods of isolated binaries before and after 10 Gyr 
of evolution by tidal friction, (a) Histogram of the assumed initial 
period distribution, (b) Histogram of the final period distribution, 
showing no dramatic change. 




Pm [d] 

Fig. 5. — The periods of the inner binaries of the simulated triples 
before and after 10 Gyr of evolution by Kozai cycles with tidal fric- 
tion, (a) Histogram of the assumed initial inner period distribution, 
(b) Histogram of the final period distribution, showing the pro- 
duction of numerous close binaries with 0.1 d < Pin, final < 10 d, 
many of which initially had much longer periods. 



This approach implies that the mass of the third star was 
correlated with the mass of the inner binary, but we do 
not believe that this correlation has any significant effect 
on our results. 

Two periods and eccentricities were picked as above. 
The smaller (larger) period was assigned to the inner 
(outer) orbit. The semi-major axes were computed from 
these masses and periods assuming non-interacting Kep- 
lerian orbits. The mutual inclination distribution of the 
tertiary is assumed to be isotropic with respect to the in- 
ner binary; thus we selected cost to be uniform in [—1,1]. 
The other angles, uj in and f2,- n , were selected uniformly 
in[0,27r]. 

After these parameters we re selected, we used the em - 
pirical stability criterion of Mardling & Aar sethl (|2001f ) 



to determine whether the system is hierarchical or if it 
will disrupt in a small number of dynamical times. If the 
semi-major axes obeyed the criterion: 

a out /a m > 2.8(1 + q out )^ ^±^tl (l - 0.3*7180°), 

(36) 

then we accepted the triple as stable and integrated its 
averaged equations of motion. Otherwise, we assumed 
it disrupted, resulting in an unbound binary and single 
star (we do not include those binaries in the following 
results). About 40% of selected triples failed to fulfill 
the condition ([36]) . A total of 7 x 10 4 stable systems 
were integrated and the results are presented here. 

5.2. Stopping conditions 

In most cases we stopped the integrations at 10 Gyr, 
roughly the main sequence lifetime of a IMq star. How- 
ever, for some systems a straight-forward integration of 
the averaged equations of motion was prohibitively ex- 
pensive. In these cases we used the following procedure 
to deduce the final state without a costly integration. 

The largest such group is triples whose Kozai cycle 
does not cause pericenter passages close enough for tidal 
dissipation to be effective. For these we integrated the 
equations until the first eccentricity maximum and com- 
pute d th e eccentricity damping timescale (Vi + V^) -1 
(eq. |A1| ) there. If it was longer than 10 Gyr, we in- 
tegrated until a second eccentricity maximum, then took 
the properties at a random time in the interval between 
the m axima, similar to the method of iTakeda &: Rasiol 
( 2005) . These systems will oscillate for their whole main 
sequence lifetimes, so choosing a random point of an os- 
cillation near the initial time is statistically indistinguish- 
able from a random time at the currently observed epoch. 

If the triple is strongly hierarchical initially, or if it is 
driven to such a state by KCTF, then the pericenter pre- 
cession due to relativity and stellar distortion dominates 
that of the third body. As shown in ij3.2[ Kozai cycles 
are suppressed in this case, so the subsequent evolution 
is very similar to that of an isolated binary. Therefore 
we stopped the integration once the Kozai timescale was 
more than 30 times the pericenter precession period, and 
evolved the system to 10 Gyr by the equations for eccen- 
tricity damping neglecting the third body (see Ej4j) . In 
real systems, the third body actually continues to have a 
modest effect which i s not modeled by the orb it-averaged 
equations of motion (jMazeh &: Shahaml [l~979). which we 
neglect. As long as this timescale criterion was satisfied, 
we also neglect the effect that nodal precession of the 
inner orbit has on tidal dissipation. 

Finally, some systems took many thousands of Kozai 
cycles to evolve significantly. We stopped individual in- 
tegrations at 4 CPU minutes, and if the system had not 
reached 10 Gyr and appeared still to be evolving, we 
re-integrated it with an artificially small viscous time so 
that the evolution would take place in ~ 100 Kozai cy- 
cles. For some individual systems we checked that the 
final state of this integration had parameters to within a 
percent of those of the final states of the original systems. 
These systems either reached the end of their allotted 
time, which was scaled down from 10 Gyr in proportion 
to the scaling of the viscous time, or stopped oscillating, 
which allowed the neglect of the third body in integrating 



10 



Fabrycky and Tremaine 




0.2 1 3 10 30 

Pbincry [d] 

Fig. 6. — The fraction of bin aries with tertiaries. T he points are 
from the observational study of Tokovinin ct al. ( 2006) — horizontal 
bars indicate the period range and vertical bars represent the er- 
ror on the fraction of tertiaries. The gray theoretical histogram is 
constructed via a linear combination of the final distributions com- 
puted in ^4] and f[5] Two free parameters were varied to achieve 
best fit with the observational results: (a) the overall fraction of 
triples relative to all systems (binaries plus triples) — 0.25 was best- 
fit — and (b) the cutoff period of the primordial distribution — 6 d 
was best-fit, corresponding to histograms in Figures |4p and |5p 
including only the two darkest shades. 



the further evolution as above. 

5.3. Numerical results 

Figure [5] shows the relation between the initial and 
final period distributions for the inner binary. Shaded 
portions of the histogram show how initial periods map 
to final periods. The main result is the strong peak in the 
distribution of periods near Pi n , final — 3d. The great 
majority of these systems have evolved onto circular in- 
ner binary orbits for which the perturbation of the third 
star no longer causes interesting effects. As an aside, we 
note that the initial period distribution of inner binaries 
(Fig. [S^) is peaked towards lower values than the ini- 
tial period distribution of isolated binaries displayed in 
Figure [4^; the latter is simply the iDuquennov fe Mavorl 
(1991) distribution. This difference is a consequence of 
the definition of an inner binary, which biases inner bina- 
ries to lower periods and outer binaries to higher periods. 

Let us pause to compare the final period distribution 
to the observed systems. In Figure|()]we plot the fraction 
of spectroscopic binaries d etermined to have a tertiary 
bv iTokovinin et all (|2006l ). We also plot a theoretical 
distribution determined by our integrations. This dis- 
tribution was constructed by taking a linear combina- 
tion of the final period histograms for isolated binaries 
(Fig. [Us) and for inner binaries of triples (Fig. Eb). Two 
free parameters were d etermined by fitting the data from 
ITokovinin et al.l (|2006f ): (a) the relative number of triple 
systems to all systems — both binaries and triples — and 
(b) the primordial cutoff to the period distributions. The 
best-fitting proportion of triple systems to all systems is 
0.25, and the best-fitting primordial cutoff period is 6 d. 

These parameters are consistent with independent esti- 
mates. The observations of additional components to bi- 
nary stars are quite incomplete, but approximately one- 



third of visual binary stars have an additional c ompo- 
nent (|Tokovinin fc Smekovl [20021 : iTokoviiiinl l2004f ) . con- 
sistent with our estimate. A complication is that the 
observed mutual inclination distribution of triple stars 
with long-period inner binaries shows moderate correla - 
tion between the inner and outer orbits (Tokovinin 1993). 
Rece nt data suggests (i) ~ 73° ±6° ([Sterzik fe Tokovinin! 
2002) , and the distribution can be reasonably represented 
by the sum of an isotropic distribution (cos i uniform in 
[—1,1]; 75% of systems) and a strongly correlated dis- 
tribution (i uniform in [0°,40°]; 25% of systems). The 
former are the systems we have simulated, and the lat- 
ter do not evolve substantially by KCTF. Therefore, due 
to this correlation, our overall tertiary fraction of 0.25 
may need to be augmented by a factor of ~ 1/3, in 
which case it is still consistent with the triple frequency 
estimated by Tokovinin. Similarly, our best-fitting pe- 
riod cutoff is consistent with the size of protostars, once 
they become dynamically stable, only contracting on the 
Kelvin-Helmholtz timescale. For example, a solar mass 
star at an age of 7 x 10 4 y r has a radius of ~ 6Rq 
(|D'Antona fc M azzitclli 1994); a binary consisting of two 
such stars in marginal contact would have a period of 
3.4 d. 

The comparison with the observed frequency of ter- 
tiaries as a function of period is an encouraging confir- 
mation of the importance of KCTF for modifying the 
period distribution of binaries. We turn now to predic- 
tions of the mutual inclination distribution, which can 
be tested by future observations. 

The relationship between initial and final inclinations 
is illustrated in Figure [7^, for systems with Pinjinai be- 
tween 3 and 10 d. During Kozai cycles, the inclination 
tends to move away from 90°. Tidal dissipation at max- 
imum eccentricity seals in these more moderate inclina- 
tions. The most interesting feature in the distribution of 
i final (Fig. U&) is the spikes that appear near i c = 39.2° 
and 140.8°. The shaded histograms show the result bro- 
ken down by initial inner binary period. Triple systems 
that start with Pi n ,initial between 3 d and 10 d, but 
do not evolve by KCTF, contribute a rather isotropic 
distribution (lightest shading). Therefore the spikes are 
fractionally much stronger if the primordial inner binary 
period distribution is cut off at a large period. This in- 
clination distribution is a distinctive feature of Kozai cy- 
cles, as predicted by the simple model of <j3](see Fig. [2})), 
and therefore would provide unambiguous observational 
evidence for KCTF. 

In a small fraction of triples, the third star is close 
enough (< 1 AU) so that its strong perturbation causes 
alignment of the inner and outer binaries even after the 
inner binary's orbit has circularized. Such systems gen- 
erally evolve to either ifi na i = 0° or ifi na l = 180° (see 
Fig. [7K, but the majority of systems that undergo this 
process have Pinjinai < 3 d). The physical mechanism 
is that the spins of the stars of the inner binary can- 
not come into alignment with the orbit normal, because 
of nodal precession due to the third body. Therefore 
tidal dissipation continues, draining energy from the in- 
ner orbit. Here Hi n is still approximately conserved (only 
a small amount of angular momentum is transfcred to 
the spins), so orbit shrink age requires tha t the m utual 
inclination decreases. See iFabrvckv et al.l f2007) for a 
detailed description of this process. If coplanar triple 
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Fig. 7. — (a) The initial and final mutual inclinations for inner 
binaries with final periods between 3 d and 10 d. Systems move 
away from 90° during the high eccentricity phase of a Kozai cy- 
cle, so tidal dissipation during that phase tends to lock in a more 
moderate mutual inclination, (b) Histogram of final mutual in- 
clination for these systems. The strong spikes correspond to the 
most probable mutual inclinations at the high eccentricity portion 
of Kozai cycles (compare to Fig. [2p). These spikes are strongest 
for binaries that have shrunk dramatically by Kozai cycles; e.g., if 
primordial binaries with Pininitial < 6 d are rare, about twice as 
many systems per 5° increment of i are expected to inhabit the 
spikes relative to other configurations. 



systems are found with i = 0°, they might also be inter- 
preted as resulting from fragmentation of a rather thin 
disk, but an observation of the purely retrograde case 
(i = 180°) would seem to require this dissipation mech- 
anism. Our simulations yield equal numbers of systems 
with i final — 0° and i final — 180°, but this is an artifact 
of our assumption of isotropic initial conditions, which is 
probably not valid for such compact systems. 

Now we step back from these small periods to survey 
the whole range of periods affected by KCTF. In Fig- 
ure [8] we plot the distribution of systems as a function 
of both cos i final and Pin, final- The spikes of Figure Wp 
are prominent at low periods, but there is another strik- 
ing feature: a deficit of long-period inner binaries with 
i w 90°. These systems have evolved to lower periods 
by KCTF. Here it is clear that KCTF can remove a sig- 
nificant fraction of near-perpendicular systems to small 
periods, over a wide range of inner binary periods, from 
~ 10 d to ~ 10 5 d. For comparison, consider an initially 
circular inner binary that undergoes eccentricity growth 
by the Kozai mechanism, then is tidally circularized from 
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Fig. 8. — Two dimensional histogram of the cosine of the final 
mutual inclination of inner and outer binaries versus the final in- 
ner binary period. Contours are spaced in increments of 50, and 
the gray-scale represents the number density with a finer grada- 
tion. Bin size is a quarter of an order of magnitude in Pi n final 
and 0.1 in cos i final, which is the same resolution as the tick 
marks. The striking paucity of systems with i final ~ 90° an d 
10 d < Pi n . fi n ai < 10 5 d is due to attrition: systems that started 
with these values have been removed by KCTF to smaller periods. 
Many of these systems inhabit the spikes of Figure [7p, which are 
clear enhancements in this plot. The region between the dashed 
lines corresponds to initial values of systems that may lose enough 
angular momentum by Kozai cycles, and enough orbital energy by 
tidal friction, to circularize to final periods less than 10 d (see text). 



Gin,max at constant orbital angular momentum. By this 
prescription, the locus of initial systems that produce bi- 
naries with Pin, final = 10 d is indicated by dashed lines 
on Figure [HI Systems between those lines may be ex- 
pected to evolve significantly by KCTF, if other sources 
of pericenter precession are negligible. The deficit of 
tri ple systems with p erpendicular orbits was identified 
by iHarringtonl (|196a ) in the first theoretical paper on 
KCTF, and this effect should be tested by determining 
the mutual inclination distribution of triple stars. 

6. PERSISTENCE OF THE FINAL MUTUAL INCLINATION 

The final product of KCTF is a triple system with a 
very large period ratio P ou t/Pin, a circular inner orbit, 
and a mutual inclination that is often near the critical 
one. Since the distinctive inclination distribution that 
results from KCTF (Fig.JTJa) may ultimately provide the 
strongest observational evidence for this process, it is im- 
portant to ask whether the inclination persists from the 
end of KCTF to the present. We point out three mecha- 
nisms for changing the mutual inclination but argue that 
each has a negligible effect. 

First, the Galactic tide may cause the outer binary 
to precess. The precession period for a binary with pe- 
riod P is roughly Pg a i/Pb, where P ga i is the vertical 
oscillation period of the star in the Galactic disk, with 
a local value of [ir/(Gp)]^ 2 « lO 8 yr[p/(O.lM p C - 3 )] 1 / 2 . 
For the precession period of the outer binary to be less 
than the age of the Galaxy, Pq w 10 Gyr, we require 
P b > Pg al /Po ~ 1 Myr. That is, the Galactic tide can 
change the orientation of the outer binary substantially 
only if a ou t > 10 4 AU. Pericenter precession of the in- 
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ner binary due to such a distant third body will likely to 
be overwhelmed by relativistic precession; equation (l34|l 
implies that TLOGR\e in =Q must be < 1 for substantial ec- 
centricity oscillations, so Kozai cycles will be suppressed 
unless ain,initicd is very large. For example, an inner bi- 
nary with an initially circular orbit can only satisfy this 
constraint if a intini ti a i > 25 AU, given a out = 10 4 AU. 
For an orbit with dinjniual = 75 AU, the pericenter dis- 
tance < 0.03 AU required for substantial tidal friction re- 
quires a maximum eccentricity in the Kozai cycle given 
by 1 — e max < 0.0004, and this can only be achieved 
if the orbits are initially almost exactly perpendicular 
(\i — 90°| < 1°). Larger or smaller values of ain,initial 
have an even stricter requirement. Therefore the num- 
ber of systems for which KCTF yields a close binary, 
after which the Galactic tide reorients the outer binary, 
is negligibly small. 

Second, individual passing stars will perturb the outer 
orbit, changing its orbital elements. The timescale for 
changing the outer binary's angular momentum is simi- 
lar to the timescale of changing its energy, which corre- 
sponds to the disruption timescale. For the typical den- 
sity (O.lM /pc 3 ) and velocity dispersion (~ 40 km/s) of 
stars in the disk, this timescale is shorter than 10 Gyr 
only for binary semi-major axes above ~ 3 x 10 4 AU. As 
above, such distant third bodies will not produce short- 
period inner binaries via KCTF with any reasonable fre- 
quency, so mutual inclination change by passing stars is 
also ineffective. 

Finally, let us specialize to a planet that has migrated 
by KCTF. Generally its orbit will not lie in the host star's 
equatorial plane, since tidal dissipation in the star will 
be too weak to align the stellar spin with the planetary 
orbit (see 97.4|) , Therefore, the planet induces a spin 
precession of its host star with a period of ~ 10 5 yrs. 
As the stellar spin precesses, the planet's orbit must also 
precess to conserve angular momentum. However, the 
magnitude of the associated change in i is < 10°, since 
the orbit of a hot Jupiter usually has more angular mo- 
mentum than the spin of its host, so this effect is unlikely 
to change the inclination distribution significantly. The 
thick appearance of the line in the final state in Figure[lJ; 
is this type of oscillation. 

7. APPLICATION TO HOT JUPITERS 

Stars are born by the fragmentation of molecular 
clouds. Planets, however, are believed to form in proto- 
stellar accretion disks, after material has settled around 
a star. In binary star systems, it might be expected 
that the angular momentum of the disk around each star 
would come into alignment with the angular momentum 
of the binary orbit. The star accretes high angular mo- 
mentum material from this orbit-aligned disk, so the spin 
of the host star of the planetary system would likely be 
in alignment with the companion orbit as well. Thus we 
might expect that both stellar spin angular momenta, 
the stellar orbital angular momentum about the system 
barycenter, and the planetary orbital angular momen- 
tum around its host star are all aligned if the binary 
semi-major axis is not too large. lHald (|1994f ) has mea- 
sured inclination to the line of sight of the spins of stars in 
binaries by comparing the rotational period of starspots 
to the wsini values of rotationally-broadened lines. If 
the stellar spins are aligned with each other, they will 



have zero difference in inclination to the line of sight. 
The converse is not true, but the prevale nce of align- 
ment can still be assessed statistically. Thus Hale (1994) 
inferred that binaries are spin-aligned for a < 30—40 AU, 
but become randomly oriented for larger orbits. There- 
fore, for nearly the entire exoplanet sample in binaries, 
it is likely that the protoplanetary disk was not aligned 
with the companion's orbit. Therefore planets arising 
from these non-aligned disks may have a high inclina- 
tion with respect to the binary companion, so KCTF 
may cause substantial evolution to the orbital distribu- 
tions of planets in binary systems. An alternative, but 
rarer, mechanism to produce non-aligned planets in bi- 
nary star systems is to form a planetary system around 
a single star, to which is later added a binary compan- 
ion through a dynamical int eraction in the birth cluster 
(jPfahl & Muterspaugh|[2006l) . 

In the past decade, about 200 giant planets have been 
discovered, of which about 20% have orbital periods short 
enough so that tides are important. The high frequency 
of close-in planets was a surprise, since planet forma- 
tion theory suggested that giant planets can only form 
beyond several AU, consistent with the current structure 
of the solar system. Therefore, a migration mechanism is 
needed to reduce the angular momentum of these plan- 
ets by a factor of ~ 10. The leading candidate is disk 
migration, in which torques between the planet and the 
remnant protoplanetary nebula tr ansfer angular momen- 
tum f rom the pla net to the gas (jGoldreich &; Tremaind 
1980; IWa7d1[l99l . There is some statistical evidence 
that planets orbiting close to one member of a wide bi- 
nary have different properties — and hence a different for- 
mation or migration history — from planets orbiting iso- 
lated stars. Among short-period radial velocity plan- 
ets (P < 100 d), the most massive planets {M Ptmin > 
2 Jupiter masses) are preferentially found in binaries 
(|Zucker fc Mazer] [2002). and these binary companions 
are cl ose enough (a out 5, 300 AU ) that KCTF may op- 
erate (jDesidera fc Barbieril 120071 ). These short-period, 
massive planets in binary stars also have lower eccen- 
tricities than short - perio d planets orbiting single stars 
(|Eggenberger et alj l200l . For these systems, KCTF 
provides an alternative migration mechanism to inter- 
actions with the protoplanetary nebula. The predictions 
of the previous sections for the period and mutual incli- 
nation distributions should still apply in the planetary 
case, and there is an additional prediction of KCTF re- 
garding the alignment of the orbital plane with the stellar 
equatorial plane (see ^7.41 below). 

7.1. Extra precession due to other planets 

In the solar system, apsidal precession is dominated 
by the gravitational perturbations from other planets; 
in multi-planet systems that are hosted by one mem- 
ber of a binary star, precession from other planets will 
compete with the pr e cession due to t he co mpanion star 
(jHolman et al.lll997l ). IWu fc Murray! (|2003l ) showed that 
the Kozai mechanism could be suppressed by small, un- 
detected masses in the HD 80606 system, since its com- 
panion star is distant so the tidal fi e ld an d precession 
rate from it are small. Ilnnanen et al.l <| 1997T ) have inves- 
tigated the instability of the solar system's giant plan- 
ets that would be caused by a companion star in an in- 
clined orbit, as a function of inclination and companion 
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mass. M arzari et al.l (|2005[ ) have investigated how a close 
(~ 50 AU) binary companion, coplanar with the planets, 
affects their scattering. The Kozai mechanism, in sys- 
tems in wh ich it operates, may als o lead to planet-planet 
scattering (jMalmberg et aLlfeOCKjl ). 

The key to why massive, short-period planets are pref- 
erentially found in binaries may lie in the competition 
between mutual planetary precession and precession due 
to the companion star, but this is a complicated process 
and the sign of the predictions is not clear. 

Take, as an example, an ensemble of planetary systems 
in which the total mass of the planets is always the same, 
but the mass assigned to individual planets varies ran- 
domly. Consider a case in which one planet is substan- 
tially more massive than the others. Initially the mutual 
perturbations of the planets suppress Kozai oscillations 
in all of them. Once moderate eccentricities and inclina- 
tions are built up in the planetary system, perhaps by 
mutual perturbations rather than the Kozai mechanism, 
the small planets become destabilized and are ejected. 
These ejected bodies no longer suppress Kozai cycles, so 
now KCTF may cause migration for the largest planet. 
Since the Kozai cycles depended on the planet ejecting its 
neighbors, this argument suggests that KCTF migration 
is more effective for more massive planets. 

As a foil, consider an ensemble of planetary systems in 
which the total mass of the planets varies randomly, but 
is always divided into some number TV of planets of equal 
mass. If the total planetary mass is small enough, Kozai 
cycles will proceed unhindered, and one or more planets 
may migrate by KCTF. The presence of the small plan- 
ets that failed to suppress Kozai cycles may not be de- 
tectable with radial velocity surveys, but the planet that 
migrated may be detected due to the larger reflex veloc- 
ity that it induces in the star. This argument suggests 
that KCTF migration is more effective for less massive 
planets. 

Since these two arguments yield opposite conclusions, 
at this time we can make no prediction regarding whether 
KCTF migration favors high or low mass giant planets. 

7.2. The period distribution of planets in binaries 

The period distribution produced by KCTF (Figure [5]) 
is quite similar to th e "pile-up" of hot Jupiters near 3 d 
(|Ford fc Rasiol l2006f ) . We have conducted the following 
statistical test to look for direct evidence of KCTF in 
the period distribution of the exoplanet sample. We 
split the radial velocity planets 3 into two groups: (1) 
si ngle planets in a s ystem with multiple stars according 
to iRaghavan et al.l (|2006h , and (2) planets orbiting sin- 
gle stars and planets in multiple-planet systems. Group 
(1) contains the systems that could exhibit Kozai cycles, 
now or in the past, since a third body is known and 
there are not other planets that could suppress Kozai 
cycles. Group (2) are the planets for which Kozai cycles 
are less likely to be present (assuming that multiple plan- 
etary systems are roughly coplanar and that they sup- 
press Kozai cycles driven by a stellar companion) . There 
may be systems with single planets a nd an undetected 
third body that causes Kozai cycles (jTakeda k. Rasiol 
2005), but this possibility simply reduces the power — 

3 Downloaded on September 15, 2006 from 

http : //vo . obspm.fr/exoplanetes/encyclo/catalog.php 
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Fig. 9. — The cumulative distributions of two populations 
of planets: those for which Kozai cycles may have been possible 
(Group 1) and those for which they are not (Group 2). Group 1 
consists of all extrasolar planets having a binary stellar companion 
but no other detected planets; Group 2 consists of all other known 
extrasolar planets. There is no statistically significant difference 
in the distributions, thus no statistical evidence for KCTF in the 
period distribution of planets. In contrast, to date the strongest 
evidence of KCTF in triple star s is the period distribution of inner 
binaries (Tokovinin et al. 2006). 



and should not bias — the statistical test. In Figure [9] we 
plot the cumulative distributions of planetary period for 
these two groups. We employ a Kolomogorov-Smirnov 
test to assess whether they are statistically indistinguish- 
able. The biggest gap between the distributions lies at 
P = 15.766 d, actually in the direction opposite from the 
expectation of KCTF; it is D n = 0.140. For n x = 27 and 
ri2 = 163, th e probability of h aving a gap at least this 
large is 72% (|Press et aL| [l992:). Therefore, there is no 
statistical evidence for KCTF in the period distr i bution 
of exoplanets. As noted earlier, iTokovinin et alJ (|2006l ) 
have employed the same test, finding different period dis- 
tributions for isolated binaries and binaries in triples with 
a probability of 0.999, which supports KCTF in triple 
star systems. 

7.3. Planetary inflation 

The energy associated with a hot Jupiter's orbit is 
about 10 times its binding energy. If the planet has 
a high eccentricity and small pericenter distance, then 
its radius will inflate from tidal heating as the or- 
bit circularizes dBodenheimer et al.l l200lHGu et"afll2003l: 
iBodenheimer et al.ll2003[ ). For instance, for HD 80606b, 
which may be under going KCTF migration currently, 
IWu fc Murravl (|2003h estimate a tidal luminosity of 
10 28 erg s _1 . For its minimum mass of 3.4 times that of 
Jupiter, energy depos ition at this rate would i nflate the 
radius by about 30% (|Bodenheimer et al.|[200ll) . The in- 
creased radius will shorten the circularization timescale 
by a factor of ~ 9, so the planetary radius should be 
included in self-consistent integrations of KCTF migra- 
tion. Since the timescale for circularizing is a very strong 
function of the radius of the planet, the maximum semi- 
major axis out to which planetary orbits are circularized 
will be diagnostic of the planetary radius at migration. 
There are three massive (> 2 Jupiter mass) exoplanets in 
binaries whose orbits are circular, two of which (Gl 86b 
and HD 195019b) have larger periods (15.8 d and 18.2 d, 



14 



Fabrycky and Tremaine 



respectively) than exoplanets with c i rcular orbits about 
single stars (IZucker fe Mazehl l2002t lEggenberger et all 
l200l . iDesidera fc Barbieril (|2007fT argue that tidal cir- 
cularization is ineffective at these large periods, but we 
suggest that this observation, if confirmed in a larger 
sample, not only supports the KCTF hypothesis but also 
may constrain the amount of inflation that is produced 
during KCTF migration. For this explanation to hold, 
the circularization timescale must remain shorter than 
the gravitational contraction timescale as the eccentric- 
ity tends towards zero. 

In the planetary case, the lower limit for the pericen- 
ter distance resulting from Kozai cycles is set by Roche 
lobe overflow. If a highly eccentric planetary orbit is 
circularized, conservation of angular momentum implies 
that the final radius of the circular orbit is twice the ini- 
tial pericenter; thus we expect that the minimum radius 
of circular orbits produced by KCTF is twice the Roche 
limit (this co nclusion is support ed by hydrodynamic sim- 
ulations; see iFaber et al J 120051 ). Indeed, the observed 
semi-major axis distribution of hot Jupiters appear s to 
cut off at twice the Roche limit (|Ford fc Rasioll2006l ). 

This simple argum ent does not acc ount for the tidal 
inflation instability of iGu et al.l (|2003l ). The planet may 
inflate enough to overflow its Roche lobe, then lose mass 
and in crease its semi-major axis. According to IGu et al.l 
(2003), starting from a large semi- major axis before in- 
flation (which is natural for migration by KCTF), the 
planet will not be disrupted but mass loss will eventu- 
ally halt, leaving the planetary orbit with a period of 1-3 
weeks and possibly with residual eccentricity. Some of 
the mass lost from the planet may be accreted onto the 
host star, leading to pollution of the star's photosphere 
with metals. Two stars that form a binary are likely to 
have the same bulk metallicity, so enhanced photospheric 
metallicity in the star hosting a hot Jupiter is evidence 
for accretion of planet material after the pre-main se- 
quence stage. Furthermore, if isotopes that are known 
to be destroyed in stars (e.g., 6 Li) are found in large 
abundance, a constrai nt on when the plane tary material 
accreted can be given (Israclian et al. 200l|). 

These effects ma y be tested for the planet HAT-P-lb 
(jBakos et a l. 2007), a transiting planet in a binary stellar 
system. It could be in the final stages of circularization 
after KCTF migration, which would lead both to some 
residual eccentricity and to an inflated radius; the latter 
is securely observed. 

7.4. Spin-orbit alignment 

In close binary star systems, including those brought 
close by KCTF, the spins of the stars are expected to 
align with the orbital angular momentum in a timescale 
short compared to the circularization timescale. How- 
ever, a star's spin probably does not align with the or- 
bital angular momentum of a hot Jupiter it hosts be- 
cause of the small mass of the planet (see Table [T]) . A 
measurement of tp, the angle between the stellar spin an- 
gular momentum and the planet orbital momentum, can 
therefore constrain the history of the system. 

It can be expected that if a planet migrated via in- 
teractions with the disk, it will remain in the same or- 
bital plane to within an angle much less than a radian. 
However, for KCTF migration, the angular momentum 
of the planet precesses about the angular momentum of 



TABLE 1 

Expected alignment properties of most triples after 
KCTF. 



alignment of: 


system type: m\-m2-m% 
star-star-star star-planet-star 


m 2 
mi 
inner 


spin / inner-orbit 
spin / inner-orbit 
-orbit / outer-orbit 


yes yes 
yes no 
no; i » 40° or 140° likely 



the inclined companion. This precession happens gen- 
erally when a body in an inclined orbit is introduced. 
Even the coplanarity of the solar system (whose net or- 
bital angular momentum is dominated by Jupiter, and 
ip = 7°) provides constraints on unseen planets in in - 
clined orbits at large distances (jGoldreich k. W ard 1972). 
As tidal dissipation takes over during KCTF migration, 
the planet's orbit eventually couples more strongly to its 
host's equatorial bulge than to the stellar companion, 
and misalignment persists as a fossil record of the earlier 
period of precession driven by the companion. There- 
fore, if planets are formed and migrate via disk torques 
to become hot Jupiters on a timescale short compared 
to the precession timescale due to the stellar compan- 
ion, then we may expect to see approximate spin-orbital 
alignment even in binary systems. On the other hand, 
if there was a period in which the planet's orbital evolu- 
tion was dominated by the binary companion, and later 
it came to be a hot Jupiter through KCTF migration, 
then non-alignment between stellar spin and planetary 
orbit will be the norm. 

We integrated a series of 1000 systems to illustrate this 
effect. Each system had identical parameters except for 
iinitiaU which ranged from 84.3° to 90° (an evenly-spaced 
grid of cosiinitiai between and 0.1). Every system had 
mi = rn-z = 1M Q , m 2 = 1O~ 3 M , and started with 
a in = 5 AU, e in = 0.1, and w m = Q in = 0°. The 
outer binary had a ou t = 500 AU and e ou t = 0. Before 
the planet has migrated at all, the Kozai timescale for 
these systems is Tmitiai = 2.36 Myr (eq. Q]). The host 
star's spin period was set to 10 d, beginning with zero 
obliquity (-0 = 0). The planet was given a viscous time 
ty = 0.01 yr (see appendix) and started with no spin. 
Once the orbit has shrunk to aj n = 0.15 AU, planetary 
orbital plane precession is dominated by the host star's 
bulge rather than the companion, so tp stops evolving; 
we call the time it takes to reach this point t a =o.i5A(7- 
The system with the smallest inclination migrated the 
slowest, with t a=0 .i5AU — 4.5 Gyr. We found that a ret- 
rograde system with initial — 90° + Si has the same value 
of ip final as a system with initial = 90° - Si, which is 
required by the symmetries of the equations of motion 
and the chosen initial conditions. Thus, the final distri- 
bution of ip from this calculation is a prediction for hot 
Jupiters that have migrated by KCTF after beginning 
with isotropic inclination relative to a companion star. 

The results are given in Figure fTUl Systems with initial 
within 3.3° of 90°, corresponding to 60% of the total, 
only undergo one Kozai cycle; i.e., &i n attains a single 
maximum, after which and are damped. Such sys- 
tems precess a fraction of a radian, therefore they form an 
orderly sequence at the smallest periods {P final < 2.5 d). 
Systems that finished at longer periods underwent many 
Kozai cycles and the orbital plane sampled many orien- 
tations. For such systems, two aspects of the evolution 
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Fig. 10. — (Top) The final period of planets after KCTF migra- 
tion versus the final stellar obli quity ip (the angle between stellar 
spin and planetary orbit) — see t|7.4l for initial conditions. These 
integrations neglected planetary inflation, which we have argued 
may affect the final periods. For final periods less than ~ 2.5 d, 
the planetary orbit reached only one eccentricity maximum, so it 
precessed less than one radian due to the stellar companion. Cer- 
tain planets in the sequence are labeled by the length of time they 
took to migrate from cti n = 5 AU to 0.15 AU, in units of the initial 
Kozai timescale r = 2.36 Myr (eq. [l]). (Bottom) Histogram of 
if). Misalignment is common, and even retrograde (if) > 90°) orbits 
are possible. This misalignment is observable in transiting planets 
through the Rossiter-McLaughlin effect and would provide strong 
evidence for KCTF migration of a hot Jupiter in a wide binary 
system. 



allow us to predict the range of ip attained. First, during 
most of these cycles, the host star maintains its space ori- 
entation, as the torque from the planetary orbit is small. 
Second, tidal dissipation is only important when ej„ is 
near its maximum, such that i is close to i m i n « 40° 
(Fig. [2]). Together, these considerations imply the sys- 
tems will produce hot Jupiters whose host stars have 



obliquities in the range i 



initial ^min 



r^j 'P r^j ^ initial ^^min- 

This range approximately describes the results for plan- 
ets with P final between 3 d and 5 d. For P final > 5 d 
the stellar spin is torqued through a large angle as the 
planet's eccentricity damps, partially realigning it with 
the planet's orbital angular momentum, sliding the range 
of final ip to lower values. (We have verified that this final 
partial realignment does not occur if the planet's mass 
is 10~ 6 Mq; its orbital angular momentum is then only 
~ 1% of the stellar spin momentum, so the planet cannot 
substantially reorient the star.) 

So far, besides the solar system, there are only four 
planetary systems for which ip has been constrained 
(see Table LI The Rossiter-McLaughlin (RM) effect 
(iRossiterl 119241: lMcLaugh!ir][l92l iGaudi fc Winnl 12007b 
constrains the sky-projection of ip (called A in the lit- 
erature) for planets that transit their host star. The 
RM effect arises because the spectral lines of the star 
are rotationally broadened, and the planet blocks por- 



TABLE 2 

Measurements of the sky-projection of the angle between 
stellar spin and planetary orbit (a) for transiting 
exoplanets. Of these, only HD 189733 is known to be a 

STELLAR BINARY. 



Exoplanet 
HD 209458b 
HD 189733b 
HD 149026b 
TrES-lb 



A 

-4.4° ± 1.4° 
-1.4° ± 1.1° 
11° ± 14° 
30° ± 21° 



Reference 
Winn et al. (20051 



Winn ct al. (2006) 
Wolf et al. (2006) 
Narita et al. (20071 



tions that are red-shifted or blue-shifted as a function of 
time when it transits the stellar disk. The derived val- 
ues of the alignment angle are effectively sky-projected 
(hence underestimated), but they are of the same mag- 
nitude as the solar system and are inconsistent with the 
distribution of ip in Figure [TUl suggesting that at most a 
modest fraction of hot Jupiters are produced by KCTF. 
A prime candidate for the measurem ent of the RM effec t 
is the transiting planet HAT-P-lb (jBakos et al.ll2007h . 
for which misalignment would be corroborating evidence 
for the KCTF migration hypothesis, made plausible by 
its inflated radius. 

Alternatively, the timing of migration can be strongly 
constrained by the p recise alignment o f the transiting 
planet HD 1 89733 b (jWinn et al.l l2006| ). According to 
iBakos et~a l. (2006) there is an M-dwarf companion at a 
projected distance of 219 AU, and preliminary measure- 
ments exclude coplanarity between the inner and outer 
orbits at the 4 — a level. If the outer orbit's misalign- 
ment is confirmed, it may be possible to constrain the 
timescale of formation and disk migration to less than 
about r sa 2 Myr, assuming the planet formed with a 
period of about 10 yr and migrated before torque from 
the companion caused misalignment. 

Further, the precession of a planet's orbital plane due 
to its host star's rotational bulge ma y be detectable 
throu gh a changing transit duration (|Miralda-Escudel 
2002). Large misalignment, which is common according 
to Figure fTTJl causes a large observable effect. 

8. DISCUSSION 

We have described the evolution and final orbital ele- 
ment distribution of binary stars resulting from the com- 
bined effects of the gravitational influence of a tertiary 
companion star and tidal friction. In particular, we 
show that Kozai cycles plus tidal friction (KCTF) can 
strongly enhance the number of binary stars with periods 
in the range 0.1 d to 10 d; briefly, the distant compan- 
ion induces strong eccentricity oscillations (Kozai cycles) 
in systems in which the inner and outer binary orbital 
planes are nearly perpendicular, and these eccentric or- 
bits are circularized near pericenter by tidal friction. 

The logarithmic period distribution of inner binaries 
produced by KCTF shows a peak near 3 d and declines 
from 3 d to 10 d. These features are consistent with 
the empirical period distrib ution of inner binaries of hi- 
erarchical systems found bv lTokovinin et ail (|2006f ). and 
we confirm these authors' interpretation that KCTF is 
the cause of this feature. Therefore KCTF appears to 
produce most close binaries. 

The dominance of KCTF for producing close binaries 
may be confirmed by testing our prediction of a modified 
mutual inclination distribution. We expect an enhance- 
ment of systems with mutual inclination close to the crit- 
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ical values of 39.2° or 140.8° if P in is in the range 3-10 d. 
Interferometric measurements have been determining the 
mutual inclination in triple stars wit h inner periods small 
enou gh for tides to be important (|Muterspaugh et all 
2006), and with improvements in limiting magnitude and 
precision they will be able to measure enough systems to 
test this theoretical prediction. Planetary inclination an- 
gles relative to a stellar companion will in some cases be 
measurable by SIM PlanetQuest. 

As the observations improve, the present theory should 
be augmented to refine the predictions. In particular, 
most of the systems in the spikes in the inclination distri- 
bution (Fig. [7b) have 2 < G ou t/L in < 10, so our assump- 
tion that the outer binary dominates the angular momen- 
tum is only marginally valid. A more accurate model for 
the known three-body gravitational dynamics may also 
have some effect on t he period d i stribu tion produced by 
KCTF. For instance, iFord et all (|2000[ ) showed that in- 
cluding octupole terms in the interaction Hamiltonian 
causes Kozai cycles to be quasi-periodic, and some peri- 
center approaches can be much closer than those com- 
puted with quadrupole terms alone (their Fig. 5). These 
closer approaches may lead to smaller final periods after 
circularization than we have predicted. Another aspect 
of the present theory that could be improved is modeling 
dynamical tides as well as, or instead of, the equilibrium 
tide. This more sophisticated tidal theory follows how 
the tidal force excites free modes of oscillation in each 
star and how these modes are damped. It may be useful 
in the present application, in which high eccentricities are 
common and the timescale of pericenter passage is com- 
parable to the periods of oscillation modes. However, the 
natural first step to refine the present analysis is to im- 
prove the model of point-mass gravitational interactions 
(which are completely understood) before attempting to 
improve the model of tidal interactions (which are quite 
uncertain). 

We have investigated the period distribution that is 
produced by KCTF alone and found that it preferentially 
produces detached binaries with Pj n ~ 2 d. However, 
over a timescale of a few billion years, magnetic winds 
may be able to extract enough angular momentum from 
such a binary s o that its com ponents come into contact 
or even merge (|Stepienl 119951) . A combination of KCTF 
and magnetic braking could therefore increase the space 
density of contact binaries. KCTF and magnetic winds 
have often been regarded as competing theories of angu- 
lar momentum loss for close binaries. However, both are 
simple additions of an extra star to well-grounded em- 
pirical relations: angular momentum loss by magnetic 
winds adds a binary companion to the spin-down pro- 
cess of single stars, and KCTF adds a third star to the 
circularization process of close binary stars. Therefore 
we believe both mechanisms must be operating at some 
level, and their relative contribution to the orbital evolu- 
tion of close binaries is an interesting subject for future 
work. 

We have considered hierarchical triple systems in the 
field of the Galaxy, where the space density of stars is 
small enough that the triples are presumably dynami- 
cally isolated (see, for instance, $5]). Of course, triples 
could exhibit KCTF while embedded in a cluster of stars, 
where the dynamics would be much richer. The outer bi- 
nary of a primordial triple would be susceptible to disrup- 



tion or reorientation by, or exchange with, passing stars 
during KCTF. Also, dynamical interactions may add a 
third body to primordial binaries (at least temporarily), 
in which case the uncorrelated orientations of inner and 
outer binaries will be above the critical inclination ~ 77% 
of the time, so Kozai oscillations will be the norm. Some 
of these effects have been c aptured by the integrations of 
iMardling fc Aarsethl (|2001h . who added tidal dissipation 
to an N-body simulation of a star cluster. A prime tar- 
get for theorists is to explain the period distribution of 
binaries in the young, nearby Hyades cluster. Its period 
distribution (for all binaries, not just inner binaries of 
hierarchical systems) shows a sec ondary peak at a few 
days that has elude d explanation dDuauennov fc Mayor! 
1 1 9 9 If) . According to lStepienl (|1995f ). angular momentum 
loss by magnetic winds cannot account for this feature — 
it tends to reduce the number of binaries at such periods, 
whereas a mechanism is needed to pile binaries up there. 
We suggest that KCTF may be that mechanism. To 
evaluate this hypothesis quantitatively, the above effects 
would need to be taken into account, as well as the clus- 
ter's young age which implies that KCTF has not yet 
progressed to completion in all cases. KCTF should first 
establish binaries with periods less than 3 d in only a few 
Kozai cycles, then slowly add longer-period binaries over 
tens to hundreds of Kozai cycles (as an example, see the 
top panel of Fig. [TU|) . This prediction may be tested by 
surveying clusters of different ages. 

Although KCTF appears to produce most close bi- 
nary stars, KCTF may produce some but not most hot 
Jupiters, for two reasons. First, systems with one planet 
orbiting a star of a binary system (i.e., those for which 
Kozai cycles are likely) have a planetary period distri- 
bution that is statistically indistinguishable from other 
planetary systems — see $7.2\ and Figure [U Second, the 
limited information on spin-orbit alignment for the four 
stars in Table [2] suggests that most host spins and plan- 
etary orbits are aligned, contrary to the prediction of 
KCTF— see flTH and Figure [TUJ For an individual bi- 
nary star system that hosts a transiting planet, however, 
large misalignment detected by the Rossiter-McLaughlin 
effect would constitute strong evidence for KCTF in that 
s ystem. 

iDesidera fc Barbieril ([2007) have shown that massive, 
short-period planets also tend to have binary companions 
within a few hundred AU, but those authors do not favor 
KCTF migration because two of those planets have peri- 
ods apparently too long for tidal dissipation to be impor- 
tant. However, during the period of extreme eccentricity 
set up by Kozai cycles, the planet's radius would expand 
due to the heat generated by tides, enhancing tidal fric- 
tion. Perhaps this effect can circularize planetary orbits 
to greater periods. To evaluate this suggestion theoreti- 
cally, planetary radius inflation must be modeled during 
calculations of KCTF. Another worthwhile theoretical 
investigation would be the the suppression of Kozai cy- 
cles by mutual planetary perturbations. 

9. CONCLUSIONS 

We conclude by reviewing the predictions of KCTF: (1) 
Period distributions for hot Jupiters and close binaries 
in systems with stellar companions have a peak at a few 
days and a minimum at longer periods (~ 10 d) for which 
tidal circularization becomes ineffective. (2) For triples 



Kozai cycles with tidal friction 
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in which the inner binary has a period between 3 and 
10 d, there is an enhancement of systems whose mutual 
inclination between inner and outer binaries is near 40° 
and 140°; for triples with longer inner periods, there is a 
paucity of systems with nearly perpendicular orbits. (3) 
Hot Jupiters resulting from KCTF migration generally 
have orbital angular momentum misaligned with the stel- 
lar spin axis by large angles, frequently even larger than 
90°, and they may have orbits circularized to greater pe- 
riods than for hot Jupiters orbiting single stars due to 
inflation of the planetary radius during migration. 

The role of KCTF in the formation of close binary stars 
and hot Jupiters can be tested by measurements of the 
following quantities in carefully defined samples of stars: 
(1) the period distribution of close binaries that have 
a tertiary component; (2) the mutual inclination distri- 
bution of triples, in order to detect the enhancement of 
systems near the critical angles (these measurements will 



normally require interferometry); (3) the angle between 
stellar spin angular momentum and planetary orbital an- 
gular momentum for planets hosted by stellar binaries, 
using the RM effect, transit timing, or determination of 
the orbit of the outer binary for transiting planets; (4) 
differential measurements of metallicity in stars in binary 
systems that host planets, to detect stellar pollution from 
overflow of the Roche lobe of the planet; (5) the distribu- 
tions of mass, eccentricity, and period for hot Jupiters. 
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APPENDIX 

TERMS FOR THE EQUATIONS OF MOTION 
Here we give the functional forms of the coefficients in the differential equations ((3][6]) . 
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Here Zj„ = 2n/Pi„ — (G(m\ + fi2)/af„) 1 ^ 2 is the mean motion. Similar expressions with subscripts 1 and 2 swapped 
hold for the second body. The parameters (X, Y, Z) form a vector in the (e^ , i\i n , hj n ) frame that gives i ts angular 
precession rate relative to the inertial frame. Their relation to the orbital elements is (|Eggleton et al.lll998[ ): 



X = i cos LOi„ + Qi n sin L0i n sin i (A6) 
Y = — i sin u>i n + cos w m sin i (A7) 
Z — Co i n + Cli n cos i . (A8) 

The dissipationless terms in these parameters can be found by computing i, Coi„ and fii„ with the Hamiltonians of £j3j 
For dissipation in mi, the tidal- friction timescale is defined in terms of the viscous timescale tvi (which we take to be 
a constant): 

tFx = ^ ( a ir) 8 < f\ a + 2^)- 2 . (A9) 

9 \RiJ (mi + m 2 )m 2 

Here, ki is the classical apsidal motion constant, a measure of quadrupolar deformabilit y which is related to the Love 
number (fet = 2ki) and the coefficient Qe given by lEggleton fc Kiseleva-Eggletonl ( 2001[): k\ = ^Qe/(^ — Qe)- W e use 



the typical value k\ = 0.014 valid for n = 3 polytropes when representing stars (Eggleton & Kisclcva-Egg letonll2001f ) 
and k\ = 0.25 valid for n — 1 polytropes when representing gas giant planets. Again, there is an analogous equation 
for fc 2 . This viscosity causes the tidal bulge to lag the instantaneous direction of the companion by a constant time. In 
the iGoldreich fc Soterl (|1966f ) theory of tides, a quality factor Q is taken to be a constant. The average fraction of the 
energy in the tide that is lost to frictional heat per radian of the orbit is and the tidal bulge lags th e instantaneous 
direct ion of the companion by an angle (2Q) _1 . When relating that angle to the time lag given by [Eggleton et al.l 
(|1998h . taking into account a stray factor of 2 mentioned in the appendix of lEggleton fc Kiseleva-Eggletonl (|2001[ ). we 
find: 
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k\ Grti\ tyi 



3(1 + 2^)2 Rf 

assuming the dissipation is dominant in body 1. Taking ty to be a constant means Q oc 



(A10) 



REFERENCES 



Bakos, G. A., et al. 2007, ApJ, 656, 552 

Bakos, G. A., Pal, A., Latham, D. W., Noyes, R. W., & Stefanik, 

R. P. 2006, ApJ, 641, L57 
Blaes, O., Lee, M. H., & Socrates, A. 2002, ApJ, 578, 775 
Bodenheimer, P., Laughlin, G., & Lin, D. N. C. 2003, ApJ, 592, 

555 

Bodenheimer, P., Lin, D. N. C, & Mardling, R. A. 2001, ApJ, 548, 
466 

D'Angelo, C, van Kerkwijk, M. H., & Rucinski, S. M. 2006, AJ, 
132, 650 

D'Antona, F., & Mazzitelli, I. 1994, ApJS, 90, 467 

De rekas, A., Kiss, L. L., & Bedding, T. R. 2007, preprint 

< |astro-ph/0703137|l 
Desidera, S., & Barbieri, M. 2007, A&A, 462, 345 
Devor, J. 2005, ApJ, 628, 411 

Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 
Eggcnberger, A., Udry, S., & Mayor, M. 2004, A&A, 417, 353 
Eggleton, P. P., Kiseleva, L. G., & Hut, P. 1998, ApJ, 499, 853 
Eggleton, P. P., & Kiseleva-Eggleton, L. 2001, ApJ, 562, 1012 
Faber, J. A., Rasio, F. A., & Willems, B. 2005, Icarus, 175, 248 
Fabrycky, D. C, Johnson, E. T., & Goodman, J. 2007, preprint 

(|astro-ph/0703418)| 
Ford, E. B., Kozinsky, B., & Rasio, F. A. 2000, ApJ, 535, 385 
Ford, E. B., & Rasio, F. A. 2006, ApJ, 638, L45 
Gaudi, B. S., Seager, S., & Mallen-Ornelas, G. 2005, ApJ, 623, 472 
Gaudi, B. S., & Winn, J. N. 2007, ApJ, 655, 550 
Goldreich, P., & Soter, S. 1966, Icarus, 5, 375 
Goldreich, P., & Tremaine, S. 1980, ApJ, 241, 425 
Goldreich, P., & Ward, W. R. 1972, PASP, 84, 737 
Gu, P.-G., Lin, D. N. C, & Bodenheimer, P. H. 2003, ApJ, 588, 

509 

Hale, A. 1994, AJ, 107, 306 
Harrington, R. 1968, AJ, 73, 190 

Holman, M., Touma, J., & Tremaine, S. 1997, Nature, 386, 254 
Holman, M. J., & Wiegert, P. A. 1999, AJ, 117, 621 
Innanen, K. A., Zheng, J. Q., Mikkola, S., & Valtonen, M. J. 1997, 
AJ, 113, 1915 

Israelian, G., Santos, N. C, Mayor, M., & Rebolo, R. 2001, Nature, 
411, 163 

Kippenhahn, R. , & Weigert, A. 1994, Stellar Structure and 

Evolution (New York: Springer- Verlag) 
Kiseleva, L. G., Eggleton, P. P., & Mikkola, S. 1998, MNRAS, 300, 

292 

Kozai, Y. 1962, AJ, 67, 591 

Malmberg, D., Davies, M. B., & Chambers, J. E. 2006, preprint 
( |astro-ph/0612041| | 



Mardling, R. A., & Aarseth, S. J. 2001, MNRAS, 321, 398 
Marzari, F., Weidenschilling, S. J., Barbieri, M., & Granata, V. 

2005, ApJ, 618, 502 
Mazeh, T., Krymolowski, Y., & Rosenfeld, G. 1997, ApJ, 477, L103 
Mazeh, T., & Shaham, J. 1979, A&A, 77, 145 
McLaughlin, D. B. 1924, ApJ, 60, 22 
Miralda-Escude, J. 2002, ApJ, 564, 1019 
Mudryk, L. R., & Wu, Y. 2006, ApJ, 639, 423 

Muterspaugh, M. W., Lane, B. F., Konacki, M., Burke, B. F., 
Colavita, M. M., Kulkarni, S. R., & Shao, M. 2006, A&A, 446, 
723 

Narita, N., et al. 2007, preprint l|astro-ph/0702707) 

Paczyhski, B., Szczygiel, D. M., Filecki, B., & Pojmahski, G. 2006, 

MNRAS, 368, 1311 
Pfahl, E. & Muterspaugh, M. 2006, ApJ, 652, 1694 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, 

B. P. 1992, Numerical recipes in FORTRAN. The art of scientific 

computing (Cambridge: Cambridge University Press, 2nd ed.) 
Pribulla, T., & Rucinski, S. M. 2006, AJ, 131, 2986 
Raghavan, D., Henry, T. J., Mason, B. D., Subasavage, J. P., Jao, 

W.-C, Beaulieu, T. D., & Hambly, N. C. 2006, ApJ, 646, 523 
Rossiter, R. A. 1924, ApJ, 60, 15 
Russell, H. N. 1928, MNRAS, 88, 641 
Stepieh, K. 1995, MNRAS, 274, 1019 
Sterne, T. E. 1939, MNRAS, 99, 451 
Sterzik, M. F., & Tokovinin, A. A. 2002, A&A, 384, 1030 
Takeda, G., & Rasio, F. A. 2005, ApJ, 627, 1001 
Tokovinin, A. A. 1993, AstL, 19, 383 

Tokovinin, A. A., & Smekhov, M. G. 2002, A&A, 382, 118 
Tokovinin, A. 2004, in Rev. Mex. AA Conf. Ser. 27, Environment 

and Evolution of Double and Multiple Stars, ed. C. Allen & C. 

Scarfe, 7 

Tokovinin, A., Thomas, S., Sterzik, M., & Udry, S. 2006, A&A, 
450, 681 

Tremaine, S., & Zakamska, N. L. 2004, in AIP Conf. Proc. 713: 
The Search for Other Worlds, ed. S. S. Holt & D. Deming, 243 
Ward, W. R. 1997, Icarus, 126, 261 
Winn, J. N., et al., 2005, ApJ, 631, 1215 
Winn, J. N., et al., 2006, ApJ, 653, L69 

Wolf, A. S., Laughlin, G., Henry, G. W., Fischer, D. A., Marcy, 
G., Butler, P., & Vogt, S. 2006, Senior Thesis, University of 
California, Santa Cruz 

Wu, Y., & Murray, N. 2003, ApJ, 589, 605 

Zucker, S., & Mazeh, T. 2002, ApJ, 568, L113 



